Kondo behavior in the asymmetric Anderson model: Analytic approach 
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The low-temperature behavior of the asymmetric single-impurity Anderson model is studied by 
diagrammatic methods resulting in analytically controllable approximations. We first discuss the 
ways one can simplify parquet equations in critical regions of singularities in the two-particle vertex. 
The scale vanishing at the critical point defines the Kondo temperature at which the electron-hole 
correlation function saturates. We show that the Kondo temperature exists at any filling of the 
impurity level. A quasiparticle resonance peak in the spectral function, however, forms only in 
almost electron-hole symmetric situations. We relate the Kondo temperature with the width of the 
resonance peak. Finally we discuss the existence of satellite Hubbard bands in the spectral function. 
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I. INTRODUCTION 

Single-impurity Anderson model (SIAM) has at- 
tracted much interest of theoretical physicists from its 
introduction in the early sixties of last century as a sim- 
ple model for formation of local magnetic momentSfi One 
of the reasons for the interest in SIAM has been its con- 
nection with the s — d exchange model and the Kondo 
effect)^ demonstrating a nontrivial strong-coupling behav- 
ior. The model won on importance after finding an exact 
solution with an algebraic Bethe ansatzi^ Revival of in- 
terests in reliable methods suitable for solving impurity 
models was brought with the concept of the dynamical 
mean- field theory of correlated lattice electrons.'* In par- 
ticular, methods addressing dynamical properties of im- 
purity models have been demanded in the mean-field de- 
scription of strongly correlated systems. Since the exact 
algebraic approaches cover only static properties, a num- 
ber of approximate analytic and numerical approaches 
have been revisited in the effort to provide an accurate 
impurity solver for strong electron correlations. 

The direct way to treat the dynamics of SIAM is to 
use the perturbation theory in the interaction strength. 
It is known that such a perturbation series converges for 
all interaction strengths and already second order gives 
rather accurate results.^^S Unrenormalized weak-coupling 
expansions break down in critical regions of any phase 
transition signaled by a singularity in a two-particle cor- 
relation function. Hence, extension of pcrturbative ex- 
pansions to the strong-coupling regime of translationally 
invariant lattice models is questionable. On the other 
hand, strong-coupling expansions based on the infinite- 
interaction model^'** fail to reproduce the Fermi liquid 
regime in the weak-coupling limit. A dynamical theory 
reliably describing the transition from the weak-coupling 
Fermi liquid to a strong-coupling solution is still missing. 

Presently, the most comprehensive quantitative ap- 
proaches producing dynamical properties of the impu- 
rity models are quantum Monte Carlo (QMC)^ and the 
numerical renormalization group (NRG)iii2i The latter 
scheme was first devised by Wilson when solving the 



s — d exchange model and later extended to SIAM in 
Ref. [HI Although designed for static thermodynamic 
properties, NRG can be used to some extent for the cal- 
culation of dynamical quantities as welli^^ii^ Since the 
quantum Monte Carlo is restricted to rather high tem- 
peratures, NRG represents now-a-days the most accurate 
quantitative low-temperature and low-energy solution of 
impurity problems for intermediate couplingsJ^ 

Unbiased numerical schemes such as NRG and QMC 
are reliable only for nonsingular quantities. A critical be- 
havior can be deduced in these schemes only from extrap- 
olations. When no simple scaling law holds we cannot 
make definite conclusions on the existence and character 
of the critical behavior in the strong-coupling regime of 
correlated electrons only from numerical data. An an- 
alytic theory is needed for this purpose. We recently 
developed an analytic approximate scheme summing se- 
lected classes of Feynman diagrams for the two-particle 
vertexii^ The approximation is justified in the critical 
region of the Bethe-Salpeter equation in the electron- 
hole scattering channel and uses a simplified version of 
the parquet equations for the electron-hole and electron- 
electron irreducible vertex functions. We demonstrated 
that in the symmetric case of SIAM, the low-energy 
Kondo scale is qualitatively correctly reproduced (uni- 
versal features thereof) in this approximation.^^ 

In this paper we discuss the strong-coupling Kondo 
behavior in SIAM in a more general context beyond the 
symmetric case. For this purpose we use simplified par- 
quet parquet equations. We introduce a general method 
of simplification of the parquet equations for impurity 
models. It is based on a partial separation of short- and 
long-range energy scales induced by a singularity in one 
of the Bethe-Salpeter equations for the two-particle ver- 
tex. In the critical region of a singularity in the vertex 
function we replace regular functions (no critical long- 
range fluctuations) by constants and keep dynamical only 
the singular ones with critical fluctuations. In this way 
we do not affect the universal features of the critical be- 
havior. Unlike the classical criticality we are, however, 
unable to uniquely separate short- and long-range flue- 
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tuations in the strong-coupling limit of correlated elec- 
trons. The universal and non-universal properties are 
mixed up there. We hence have to decouple the short- 
and long-range fluctuations with an additive approxima- 
tion to make the ansatz of separation of scales consistent. 
We can do it in several ways. Here we discuss two of 
them. In the first one, introduced in Ref. Ilisl . we demand 
fulfilling of the parquet equations for regular vertices av- 
eraged over short-range fluctuations to obtain constants. 
This is a mean- field- type approximation that can be ap- 
plied in the whole region of the input parameters also be- 
yond the critical region of the Bethe-Salpeter equation. 
It is a fully consistent approximation in the Fermi-liquid 
regime. Another, more accurate way of decoupling short- 
and long-range fluctuations introduced in this paper is to 
use systematically the dominant small energy scale shap- 
ing the singularity in the vertex function to reduce the 
frequency dependence of regular functions. The latter 
approach is more precise as concerns the critical behav- 
ior but less flexible as concerns applicability outside the 
critical region. 

Studying the Kondo behavior in SIAM we distinguish 
its two different aspects. First, we introduce a small 
dimensionless scale connected with a singularity in the 
two-particle vertex that we call the Kondo scale. From 
this scale we derive the Kondo temperature as a cutoff 
screening the actual singularity that cannot be reached 
in impurity models. Next we analyze the behavior of 
one-electron functions in the Kondo limit, i. e., with a 
quasi-divergent two-particle vertex. We demonstrate the 
existence of a narrow Kondo resonance peak near the 
Fermi energy for small deviations from the symmetric 
situation. We calculate the position of the Kondo peak 
and estimate its width. We also discuss the existence of 
the satellite Hubbard bands in the strong-coupling limit. 

It is more elaborate to study asymmetric situations 
in impurity models than the symmetric case. We must 
introduce an effective chemical potential controlling the 
occupation of the impurity level. To guarantee consis- 
tency of approximations, the chemical potential must be 
determined self-consistently from the actual filling of the 
impurity level. It means that independently of the form 
of one-electron propagators we use in the approximate 
treatment we have to calculate impurity filling from the 
full spectral function with the self-energy obtained in 
the parquet approximation. We hence must achieve a 
static one-particle self-consistency. Only in this way we 
avoid spurious (Hartree-type) transitions in the strong- 
coupling regime. 

The paper is organized as follows. In Sec. |TT] we in- 
troduce the parquet equations for impurity models. The 
simplifications of the parquet equations in the critical re- 
gion of a singularity in the Bethe-Salpeter equations are 
discussed in Sec. lIIII Two-particle functions in the Kondo 
regime are calculated in Sec. lIVI and the one-particle func- 
tions in Sec. El 



II. PARQUET EQUATIONS FOR IMPURITY 
MODELS 

We start our investigation with the Hamiltonian of 
SIAM that reads 

kcr cr 

+ E (^k 4ck. + V:clj„) + Unpf . (1) 

kcr 

We denoted = d^^da- When calculating the grand 
potential and thermodynamic properties of the impu- 
rity site we can explicitly integrate over the degrees of 
freedom of the delocalized electrons. To this purpose 
we standardly replace the local part of the propaga- 
tor of the conduction electrons by a constant A(e) = 
TT^j. |VkP^(e — e(k)) = A the value of which we set as 
the energy unit. With this simplification we expand all 
dynamical and thermodynamic quantities in powers of 
the interaction strength U and rearrange appropriately 
the resulting series. 

The principal idea of the parquet approach is to de- 
rive the self-energy and other thermodynamic and spec- 
tral properties via a two-particle vertex F. This ver- 
tex is standardly represented via Bethe-Salpeter equa- 
tions. The Bethe-Salpeter equations express the full two- 
particle vertex by means of irreducible ones that can be 
generically represented asi^ 

r = A"-f [A^GGla^F (2a) 

where A" is the irreducible vertex in channel a, the 
brackets [...]« stand for an appropriate interconnection 
of internal variables in the chosen scattering channel a 
and * represents summation over intermediate states cre- 
ated during the scattering process. We use in our con- 
sideration only singlet electron-hole and electron-electron 
scattering channels. 

The consistency of approaches with more two-particle 
vertices is guaranteed by the parquet equation that reads 

F A'^'* + A'^'^ - X . (2b) 

We introduced a vertex X irreducible in both channels 
that in the parquet approximation is replaced by the 
bare interaction U. The parquet equation (|2b[) says 
that the full two-particle vertex is a sum of the irre- 
ducible vertices from the electron-hole channel (A^'') and 
from the electron-electron channel (A'^'^) from which the 
completely irreducible vertex (bare interaction) is sub- 
tracted. Equation (j2b[) is used in each Bethe-Salpeter 
equation Pa|) to exclude the full vertex F and to close 
the approximation for the irreducible vertices A'^'' and 
A^"^ . In the impurity case the dynamical variables in the 
vertex functions are Matsubara frequencies. Hence the 
explicit form of the parquet equations with the electron- 
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hole and electron-electron singlet channels read 

- ^ Afl{iuJn, ; il'„i)Gl{iU!n")Gl{iuJn"+m) 

' n" 

-X|j,(it^„'>,iw„'ii/„i)] . (3a) 

and 

All (ia;„ , zw„' ; i ) = I|x (ia;„ , iw„/ ; w^) 

n" 

X G|(i[x)„+„'+m_„") [A||(ia;„", ZWn'; Wm+n-n") 
+A||ilx'n" , ZCJn' ; i^'m+n— n" ) 

-iMi^n" ,i^n']iVr,i+n-n")] ■ (3b) 

These equations form a set of coupled nonlinear inte- 
gral equations that cannot be solved by available analytic 
means. Their numerical solution with 2" = J7 is available 
only for high temperatures^ in the region irrelevant for 
quantum effects induced by strong electron correlations. 
Since we do not know the exact form of the completely 
irreducible vertex T we have to resort in the parquet ap- 
proach to approximations anyway. 



III. TWO-PARTICLE CRITICALITY AND 
SIMPLIFIED PARQUET EQUATIONS 

For most of physically interesting situations it is not 
necessary to know the detailed form of the completely 
irreducible vertex I. If it is a regular function it is rea- 
sonable to replace it with the bare interaction, which 
is standardly done in the parquet approximation. It is 
then futile to try to solve the parquet equations with a 
tremendous effort in regions where no significant devia- 
tions from Fermi liquid or gas are present. The principal 
area of application of the parquet approach should be 
in the low-temperature strong-coupling region where the 
standard weak-coupling schemes fail. 

The critical region of a singularity in one of the Bethe- 
Salpeter equations is the situation where the difference 
between the parquet theory and other weak-coupling 
schemes becomes significant. Parquet equations due to 
their nonlinear character introduce a two-particle self- 
consistency enabling to handle properly the critical be- 
havior. In particular, the two-particle self-consistency ef- 
fectively suppresses spurious non-physical singularities of 
one-particle approximations. That is why we try to solve 
the parquet equations, at least in an approximate way, 
in the critical region of a singularity in Bethe-Salpeter 
equations. Fortunately, singularities with divergent ver- 
tex functions offer a natural way of a simplification of the 
parquet equations. This simplification is a separation of 



long-range fluctuations of divergent quantities from the 
short-range ones of regular functions. In the critical re- 
gion of a singular vertex we then treat finite differences 
in regular functions as negligible with respect to macro- 
scopic changes shaping the singular behavior. 

To make the simplification of the parquet equations 
effective we must rearrange the weak-coupling expansion 
in such a way that we first reach the critical region of 
a singularity in Bethe-Salpeter equations. In SIAM it 
is achieved in the random-phase approximation (RPA), 
that is, the Bethe-Salpeter equation in the electron-hole 
channel with the bare interaction and the Hartree one- 
electron propagators. The two-particle vertex in Mat- 
subara frequencies then reads 



-■B.PA 



U 



1 + UXeh{iVm) 



(4) 



where we denoted Xeh{ivm) = 
l3~'^^^G](iLjJn)Gi{iix)n+m) the dynamical electron- 
hole bubble. It is easy to evaluate the electron- 
hole bubble with the Hartree propagators 
G{x + iy) = l/(a;-7I-fzsgn(y)(A-K I2/I)) to be 



Xeh{z) 



r ^/(.) 

J-00 



1 



2: — /i + X + «A 
A 



X — ji — z + i A 



[(a;-7i)2 + A2]2 



(5) 



where Jl = Ed + Un is an effective chemical potential 
measuring deviations from the symmetric case JI = 
and f{x) = 1/(1 + exp{/3x}) is the Fermi function. We 
used < rt < 1, the actual occupation per spin of the 
impurity level. Since the static value Xeh = Xeft,(0) < 0, 
we choose a finite interaction strength U and a very low 
(zero) temperature T so that we are very close to the sin- 
gularity in vertex F^^"^ with < 1 + U Xeh ^ 1- We now 
rearrange the expansion beyond RPA in such way that 
we remain all the time in the critical region of this singu- 
larity. We then replace the interaction strength U as an 
expansion parameter by a new small scale a — \ -\- Uxeh- 
We keep this scale fixed during summation of diagrams 
and use it to separate large (relevant) from small (irrele- 
vant) quantities. We will discuss ways how we can use the 
existence of a small (large) scale to simplify the parquet 
equations. In fact we will see that it is not the scale a it- 
self but rather 1/| lna| that is taken as a small parameter 
in the expansion beyond RPA in the critical region of the 
electron-hole vertex. By the parquet equations we effec- 
tively replace the direct summation of diagrams beyond 
RPA and their renormalization in the critical region of 
the singularity in the Bethe-Salpeter equation from the 
electron-hole channel. In analogy with classical critical- 
ity, RPA corresponds to a static mean-field approxima- 
tion and the parquet approach to renormalization-group 
dynamical corrections to the mean-field critical behavior. 
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A. Finite temperatures: mean-field approximation 

By solving the parquet equations in the critical re- 
gion of a Bethe-Salpeter equation we assume that no 
other critical point emerges beyond that found in the one- 
channel (mean-field) approximation. The singularity in 
SIAM emerges in the electron-hole scattering channel and 
hence the irreducible vertex from the electron-electron 
channel A||(ia;„, zct;„', w„j) is singular in the variable 
ii^m = 0. As no other critical scattering appears in the 
model, the irreducible vertex from the electron-hole chan- 
nel A|j'(iw„, iwn', Wm) remains regular. We now neglect 
all finite differences in regular functions and replace A'^'' 
with a constant, an effective interaction U. Inserting this 
ansatz in Eq. (j3ap we obtain 



A'i'l{iujn,iuJn';ivm) = U - - — =■ 



— 2 



1 + UXehiiVm) 



(6a) 



Being in the critical region where the denominator on the 
right-hand side of Eq. ^ is very small, we can neglect 
the contribution to A'^'^ from the bare interaction and 
take into consideration only the low-frequency singular 
part of this vertex 



A^™^(w„i) 



— 2 

1 + UXeh{il^m) 



(6b) 



Neglecting finite differences due to the frequency depen- 
dence of the regular vertex K'^^ we have won a simplified 
frequency dependence of the singular vertex in form of 
the RPA vertex r-^^^ from Eq. (g]). This simple form 
enables us to control analytically the critical behavior of 
the singular vertex A*^"^. 

We now insert the derived singular part of the ir- 
reducible vertex from the electron-electron channel into 
Eq. ((3b|) and obtain 



m+n' — n" 



X G't(iw„")Gx(iw„+„'+m-n") 

X [c7-f A^™»(zz.™+„_„.0 



This equation is evidently inconsistent, since its right- 
hand side is frequency dependent. We neglected finite 
differences in the regular electron-hole irreducible ver- 
tex on the left-hand side of Eq. ([7]). We must do the 
same on the right-hand side as well. It contains, how- 
ever, singular functions. We hence cannot demand equal- 
ity of both sides for a specifically chosen Matsubara fre- 
quency. To make our simplifying ansatz consistent we 
have to first regularize the right-hand side of Eq. ([7]) 
and grant equality only in a mean. To achieve consis- 
tency we multiply both sides with G\(iLiJri)Gi(iLOm-n) 
and G'|(iix'„')G|(jix'm_„') and sum over fcrmionic Mat- 
subara frequencies zw„ and iwn' ■ We further put m = 



to keep the effective interaction real. Since the low- 
frequency contributions are the most relevant ones, we 
obtain 



UXee = U Xe 



Xee + (GtCxLTi) ' 



(8) 



where Xee — P "'^ X^n '-'T(*^n)^i(*^-") static 
electron-electron bubble. Further on we used a notation 



L|x(iw„) = ^^G|(icj„OG'x(iw-„')A' 



(9a) 



{G^GiX) = i ^ G^{^u;r.)Gl{^u;.n)X{lu;n) ■ (9b) 



Equations ([6]) and ([8])- ([9]) form a closed set of rela- 
tions determining the static effective interaction U and 
the dynamical vertex K'^™^{ivm) as functionals of the 
one-particle propagators G^ and the bare interaction 
U . They were derived and justified in the critical re- 
gion of a singularity in the Bethe-Salpeter equation from 
the electron-hole channel. That is, the denominator on 
the r.h.s. of Eq. ([6]) approaches zero. The form of the 
defining equations for A'*'"'^(wm) and U enables their 
application in the whole region of the input parameters. 
The approximation remains consistent everywhere in the 
Fermi-liquid regime. 



B. Zero temperature: low-frequency approximation 

Separation of large scales (potentially divergent func- 
tions) and small scales (regular functions) is incomplete 
in SIAM. To stay close to the critical point of the singu- 
larity in A^'"^ we need to know the value of the effective 
interaction U . We do not need to know, however, the 
detailed frequency dependence of vertex A^^ but only its 
action in Bethe-Salpeter equation ([5a|) . With an effective 
static interaction we make a replacement 



^] ■ C^) (3 ^^^^\{^UJn,^^n'']Wrn)Gy{^i^n'')Gl{iUJn''+rn) 



X Kf^{iuJn" ,iuJn';iVm) 

UXeh{Wm)A'\\{iuJn,ii^n']iVyn) 

in the exact parquet equation. This replacement is qual- 
itatively correct if we treat integrals (sums over Matsub- 
ara frequencies) of singular functions as regular ones and 
as far as the effective interaction U remains finite. We 
will see that in SIAM it is always the case even in the 
limit U oo. 

Although the critical behavior of the singular vertex 
^smg j^Qgg t^qi depend on the value of the effective inter- 
action J7, we need its value to determine the asymptotic 



5 



behavior when the critical point is approached. We need 
to determine how fast the small scale 



IV. KONDO SCALE AND KONDO 
TEMPERATURE 



a = l + Ux. 



eh 







(10) 



approaches zero as a function of the input parameters, in 
particular, of the bare interaction strength U . It is evi- 
dent that this result depends on the function dependence 

uiu). 

The averaging of Eq. ([7]) we used in the preceding 
subsection to reach consistency in the determination of 
the effective interaction is only one, mean-field way. We 
can proceed in a more accurate way if we fully utilize 
the dominance of the low-frequency asymptotics of the 
singular vertex. We can approximate the singular vertex 
in the critical region with a <C 1 as 



A'"'"f(w+) 



U 



(11) 



where we denoted x' — i'!^dxehioJ+) /duj\uj=o and (jJ+ — 
lli + iO+. It is clear that only frequencies of order an/Ux' 
are relevant in integrals with this singular function. It 
means that only frequencies of order a/Ux' are relevant 
in regular functions when integrated with the singular 
vertex. We hence replace all values of the redundant 
frequencies in regular functions by zero. This is possi- 
ble only at zero temperature with a continuous distribu- 
tion of Matsubara frequencies. The nonsingular electron- 
hole irreducible vertex A'^''(iu;„, ia;„'; w„i) = U. Such 
a static replacement cannot be done, however, in sin- 
gular functions. Before we can set redundant frequen- 
cies to zero in Eq. ([7]) we multiply both sides either by 
G'|(«w„)Gx(iWm_„) or by Gi^{iuJn')G i{iuJm,-n') and inte- 
grate over iu)n or iwn' , respectively. The output is a 
" sufRciently" regular function with maximally logarith- 
mically divergent terms proportional to |lna| where we 
already can put all external frequencies to zero. We then 
obtain another equation for the effective interaction 



UXee = U Xe 



7^2 



l + L 



n 



(12) 



where we denoted L-^ = L||(0). This approximation is 
applicable only at zero temperature with the sums over 
Matsubara frequencies replaced by integrals. In this ap- 
proximation, however, the singularity in the two-particle 
vertex is treated more accurately, since unlike the pre- 
ceding construction we treated here logarithmically di- 
vergent scales of order |lna| more sensitively. Logarith- 
mic terms may go lost or be washed out by additional 
summation over fermionic Matsubara frequencies used 
in the preceding subsection. It means that our precision 
is determined by the largest "regular" scale being now of 
order | lna|. There is no such a small parameter in the 
construction with Eq. ([5]). 



We now proceed with solving the derived simpli- 
fied parquet equations for the singular vertex A*™f (ijj+) 
and the effective interaction U. We resort to the spin- 
symmetric case G| = G^ and analyze the solution in the 
critical region of the singularity with a — 1 + Uxeh — 
1-U /_!^ (iw/(w)S[G+(w)2]/7r ^ 0. We used an abbre- 
viation G±(ll}) = G{u! ± iO"*"). The new small dimension- 
less scale a emerging in the critical region of the two- 
particle vertex is dominant for the determination of the 
behavior of SIAM in the strong-coupling regime. We call 
it a Kondo scale, since, as we will demonstrate below, 
it gives origin to and controls the Kondo strong-coupling 
asymptotics known from the exact Bethe-ansatz solution. 

We first evaluate the dominant contributions from the 
singular vertex A*™s(ci;_(-) to the integral L^i in the limit 
a — > 0. We simultaneously use both reduction schemes, 
Eq. ([8]) and Eq. (fT2|) . The former for non-zero temper- 
atures and the latter at zero temperature. Using the 
low- frequency asymptotics, Eq. (fTTj) . we obtain at finite 
temperatures 



G{z)G{-z) 



dx 



tanh(6.T/2) a? + x"^ 



(13a) 



where we denoted a dimensionless inverse temperature 
6 — 7r/3/[/x', effective bandwidth ^ — Ux' D/-K, and x' = 
—4 (ia;/(Aa;-t-/l)a;(x^-f 1)"'^. At low temperatures, at 
which the integral becomes singular and the asymptotic 
representation (I13ap holds, we further obtain 



G{z)G{-z) 



2X' 



In 



1 + 62a2 



4 / 1 

— arctan — 
ba \ ba 



(13b) 



The leading singular contribution to the integral L at 
zero temperature comes out rather simply 



G{z)G{~z) 



(13c) 



where po = A/7r(A2 + is the density of states at the 
Fermi energy. 

Parquet equations (O, ([5]) do not allow for a solu- 
tion with a critical point, a = 0, at a finite temperature, 
that is with with ba — 0. The critical point can only 
be asymptotically approached when the temperature ap- 
proaches zero. Actually, the critical behavior saturates 
when the product ba becomes of order unity. This condi- 
tion in the asymptotic limit 6 — s- cx) defines a temperature 
at which the static two-particle vertex function asymp- 
totically approaches its maximal value. In the Kondo 
regime, below this saturation temperature, called Kondo 
temperature, physical quantities are saturated and the 
Kondo temperature overtakes the role of the tempera- 
ture scale. Using the low-temperature asymptotics from 
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Eq. (jl3bp together with the condition 6a = 1 in Eq. 
we obtain an expression for the Kondo temperature 



X' 



X cxp 



\Xeh\ {[GG)ie) 



(14a) 



where we denoted {{GCy^^f 

poo 

-co 



-TT 



^ (ia;/(a::)Q [(G+(a::)G_(— x))"]. Kondo regime 
sets in only asymptotically in the limit J7|Xe?i(0)| — > oo. 
The genuine strong-coupling regime in SIAM sets in 
already for ?7|xeh(0)| > 1. This regime is not reached 
by simpler weak-coupling approximations such as RPA 
or FLEX, where always U\xehiO)\ < 1- Notice that the 
Kondo temperature cannot be defined unambiguously 
and Eq. (jl4ap is valid as order of magnitude only. 

Using Eq. we cannot determine the Kondo tem- 
perature directly, since it is applicable only at zero tem- 
perature. We can nevertheless determine the Kondo scale 
a as a function of the bare interaction strength. Insert- 
ing the asymptotic result from Eq. I\l'dc\i into Eq. we 
obtain an equation for | lna|. From it we find a solution 
for the Kondo scale 



exp 



\Xeh\ 



mx. 



eh\ 



(14b) 



Expressions (|14ap and (|14bp for the Kondo tempera- 
ture and the Kondo scale were derived for an arbitrary 
form of the one-electron propagators. The only assump- 
tion was that the electron-hole bubble Xeh{z) is analytic 
at z = 0. We argued in Ref. [H that the Hartree one- 
electron propagators deliver the best results in SIAM on 
both small and large energy scales. Any dynamical one- 
electron self-consistency negatively interferes in the con- 
trol of the two-particle singularity in the Bethe-Salpeter 
equation won by the parquet approximation. In particu- 
lar, the electron-hole bubble Xeh with [/ — > oo when 
calculated with the fully renormalized one-electron prop- 
agators, which forces the effective interaction to diverge. 
Our simplifying scheme was justified for finite effective 
interaction U < oo. The parquet approximation for the 
vertex functions with T — U corresponds to a Hartree 
approximation on the one-particle level. Going beyond 
Hartree propagators in the parquet scheme demands to 
renormalize adequately the completely irreducible ver- 
tex T to keep the approximation consistent. The only 
peremptory one-particle self-consistency is the static one, 
where the effective chemical potential /I is calculated from 
the fully renormalized one-electron propagator. 

We hence use the Hartree propagators in our parquet 
approximation to evaluate explicitly the two-particle 
functions in Eqs. (jl4p . In the low-temperature limit the 
dominant contribution to the Kondo temperature comes 
form the zero-temperature values of the functions on the 
right-hand side of Eq. (|14ap . At zero temperature the 



two-particle bubbles read 
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Further on we obtain 
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The effective interaction in the Kondo regime reads U — 
^(l-a)(/l2 + A2)/A. 

With the above representation we obtain an explicit 
dependence of the Kondo temperature on the interaction 
strength U and the effective chemical potential Jl. We 
explicitly mention only two important limits. First, in 
the symmetric case Ji — Q we have 



ksTK 



A cxp 



(17a) 



The Kondo temperature is determined by the Kondo di- 
mensionless scale derived in Ref. 

The other limit with a simple representation of the 
Kondo temperature is an almost empty impurity level 
where /I S> A. In this case we obtain 



kb-Lk = exp 



27rA 



3/i \TTI1 



UA 

2 
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(17b) 



where 11 = Ed 



Un. 



For /i 3> A we explicitly have 
ji = Ed{l - ^1 + 4C/A/7r£'2)/2 and n ~ A/nJl. From 
which we get JI ^ \Eci\ and E"^ ~ UA. We see that 
the strong-coupling regime UA > ttE^ sets in also for an 
almost empty impurity level, but the Kondo temperature 
is no longer exponentially small. 

The Kondo scale calculated from Eq. (jl4bp with 
Hartree propagators reads 



7r/2 + arctan( ^^~^ ) 



a = exp 



271A 



-A-A2-7I2 

TT 



(18a) 



In the symmetric case and for an almost empty impurity 
level it reduces to 



a = exp{-(C/po - 1)} 



(18b) 
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and 



exp 



2\Eh\A 



{UA - 



(18c) 



respectively. Using further the solution for in the 
limit n we obtain an asymptotic dependence a ~ 

exp{-^|£;d|} or a - exp|-a'y^| for U oo. The 

Kondo regime sets in for almost empty impurity level if 
Un^ > TT^A. 



Both expressions (|17ap and (jl8bp for the Kondo tem- 
perature in the symmetric case qualitatively agree (in 
universal features) with the exact formula ksTx = 
y/UA/2exp{T:EdiEd + U)/2AU} known from the s - 
d exchange Hamiltonian asymptotically reached in the 
limit t/ — > oo of SIAM.^* The Kondo temperature for 
an almost empty impurity level, Eq. (jl7b[) . deviates in 
the dependence on Ed from the exact result. The two 
asymptotic results were, however, derived under different 
conditions. The former representation was obtained for 
a fixed /I S> A while the latter for a fixed Ed constrained 
to Ed < -A a.nd Ed + U > AM The Kondo scale 
from Eq. (|18c|) improves upon this imperfection in the 
asymmetric situation of the approximate scheme based 
on Eq. ([8]) . It qualitatively well reproduces the exact 
result with the exponent linearly proportional to the po- 
sition of the impurity level Ed- As expected, the zero- 
temperature, low-frequency approximation on the effec- 
tive interaction U, Eq. ((T2)) . better reproduces the de- 
pendence of the Kondo scale on the interaction strength 
and in particular on the position of the impurity level 
and the impurity-level occupation. 



V. SELF-ENERGY: KONDO RESONANCE AND 
HUBBARD SATELLITE BANDS 

We defined the Kondo temperature and the Kondo 
scale from a critical behavior of the two-particle vertex 
near a singularity in the electron-hole correlation func- 
tion. The one-electron propagators were used as input 
parameters for the two-particle vertices. In this situa- 
tion with one-electron propagators chosen independently 
of the two-particle ones we cannot interpret the two- 
particle functions we use as the actual physical quan- 
tities they diagrammatically stand for. We have to con- 
form with conservation laws in consistent theories. To 
make the approximation thermodynamically consistent, 
we have to introduce the self-energy as a functional of 
the two-particle vertex. The physical quantities in con- 
serving theories are then derived from the self-energy via 
appropriate (functional) derivatives. 

We use the Schwinger-Dyson equation to relate the 
self-energy with the two-particle vertexii^ In our case of 
the simplified parquet equations it reduces to 
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FIG. 1: Density of states for [7 = 4 and various temperatures 
demonstrating formation of the Kondo resonance at the Fermi 
energy for the symmetric case. The zero-temperature value 
of DOS in the half-filled case due to the Friedel sum rule is 
l/vr. 




1 + UXeh{iOJn' - iuJ„) 



(19) 



FIG. 2; Temperature dependence of the density of states for 
a partial occupation of the impurity level n = 0.4 and U — 4. 



In the spin symmetric case and in the Kondo regime (crit- 
ical region of the singularity in the two-particle vertex) 
we again use the low-frequency expansion of the denomi- 
nator in the analytically continued version of Eq. (fT9|) . 
We further employ the electron-hole symmetry of the 
Hartree propagators G(ia;„ — /l) = — G(— iw„-|-/2) to guar- 
antee analytic properties also in approximate numerical 
evaluations of the self-energy. After analytic continua- 
tion we explicitly obtain for the real and imaginary parts 



8 



of the self-energy 

my)- 



dx 



1 



27r(l - a)x' J-oo {x + y - -pY + 1 x"^ + a"^ 
'x{x + y-'p) 



tanh(/5a;/2 
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Stanh {f3{x + y)/2) 
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(20a) 



27r(l - a)x' J-oo [x + y-JlY + l x^ 
1 



tanh(/3a;/2) 



tanh {P{x + y)/2) 



(20b) 

where we set A = 1 and denoted a — a7r|xeh(0)|/(l — 
a)x' ■ Representation is meaningful for < a < 1, 
but it is justified in the asymptotic limit /3 — s- cx) and 
Pa ^ 1. The physics in the Kondo regime of SIAM, 
including the effective chemical potential /I, in this ap- 
proximation is determined from self-energy (|20p and its 
dependence on external sources entering the theory via 
the Hartree propagators. For a numerical solution, not 
deep in the Kondo regime, it is straightforward to iter- 
ate Eq. ^ with either Eq. ^ or Eq. ^ and use their 
solution in Eq. (jl9p to determine dynamical and thermo- 
dynamic properties of SIAM. 

Figure [T] shows temperature dependence of the spec- 
tral density for J7 = 4 and half filling. When the tempera- 
ture decreases, quasiparticle states are first expelled from 
the Fermi energy to satellite peaks. But at temperatures 
of order of the Kondo one, a new resonance around the 
Fermi energy starts to develop so that a metal-insulator 
transition is screened. This scenario takes place at ar- 
bitrary filling, but the exponentially narrow Kondo res- 
onance realizes practically only for fillings close to the 
symmetric case, Fig. [5J Farther apart from the half- filled 
case the Kondo temperature increases, the quasiparti- 
cle peak broadens and gradually merges with the lower 
satellite band. Fig. [31 Moreover, the weight of the Kondo 
peak also decreases with doping the symmetric case with 
either holes or particles. The numerical results are in 
a good agreement with the behavior obtained by NRG 
in Ref. [l^. It is important that the simplified parquet 
approximation used here predicts qualitatively correctly 
formation of a Kondo peak near the Fermi energy also 
in slightly asymmetric situations where the Friedel sum 
rule (obeyed by this approximation) pinning the zero- 
temperature density of states at the Fermi energy no 
longer holds. The Kondo peak in asymmetric fillings is 
hence a consequence of genuine strong electron correla- 
tions. Only theories producing a Kondo resonance peak 
in the asymmetric case can be relied upon in the strong- 
coupling regime. 

We used Eq. ([5]) for the determination of the effec- 
tive interaction U in the calculations of the self-energy 
and the spectral function. The differences in the results 
obtained with the other approximate decoupling of the 
parquet equations in Eq. (jl2p are not dramatic. The re- 
sults are qualitatively the same with minor quantitative 




FIG. 3: Dependence of the Kondo resonant peak on filling 
of the impurity level for U = 4. Merging of the quasipar- 
ticle peak with the lower Hubbard band for lower fillings is 
apparent. 



U = 8 



O 

Q 




FIG. 4: Spectral function calculated with the effective interac- 
tion calculated from Eq. ^ and Eq. JTS]) for U = 8. The two 
approximate schemes do not significantly differ at any filling 
except for a narrow asymptotic critical region of asymmetric 
models. 



deviations as demonstrated in Fig. |31 In the symmet- 
ric situation the distance between the satellite peaks is 
smaller and the central peak is slightly broader when 
Eq. in]) is used. 



Critical behavior of the two-particle vertex in the 
Kondo regime (a ^ 1) allows one to assess the self- 
energy and the resulting spectral function analytically. 
Using the fact that a ^ in Eqs. (PO)) we obtain for the 
leading asymptotic behavior of the self-energy at zero 



9 



U = 8 




FIG. 5: Real part of the self energy from the numerical solu- 
tion (N) and the asymptotic form from Eq. (|21|l (A) at zero 
temperature and two fillings. The critical behavior is more 
pronounced in the asymptotic solution 



temperature 



3fiS(tj+) = 



C/ln 
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— 7n — 



U arctan ( 
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C/ln 



2C/7r2po 
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Representation ([21]) is quite general and holds in the 
Kondo regime for any form of the one-electron propa- 
gator Q used in the parquet approximation. That is, we 
can analyze with this representation any dispersion rela- 
tion or approximate form of the one-electron propagators 
used in the Schwinger-Dyson equation (|19p . Most impor- 
tantly, we can assess the width of the Kondo resonance 
near the Fermi energy and decide about the existence of 
the satellite upper and lower Hubbard bands. In Fig. [5] 
we compared the real-part of the self-energy calculated 
with the full representation, Eq. pp]) , and the asymptotic 
form, Eq. (HI]). We can see that there are no relevant 
(qualitative) differences in the results and the agreement 
is almost perfect in the critical region and for small en- 
ergies near the Fermi energy. The asymptotic formula 
should hold in the critical region with a <^ 1 and hence 
the deeper in the critical (Kondo) region we are the bet- 
ter the agreement becomes. 

We already proved in Ref. [1^ that one-particle self- 
consistent approximations with fully renormalized one- 
electron propagators smear out the satellite Hubbard 



bands. On the other hand, numerical calculations clearly 
demonstrate that the approximation with the Hartree 
one-electron propagators leads to satellite bands. With 
the asymptotic formulas (j2ip we can even prove it and 
determine the positions of these bands. 

We use Hartree propagators in Eqs. (pij) . The effec- 
tive chemical potential JI is determined from the actual 
occupation of the impurity level. We first estimate the 
position of the Kondo resonance peak that is, maximum 
of the spectral function very close to the Fermi energy, 
within a frequency region of order of the Kondo scale. 
We determine the center of the Kondo peak from vanish- 
ing of the imaginary part of the one-electron propagator. 
Maximum of the spectral function at the Fermi energy 
is guaranteed by the Friedel sum rule only at half fill- 
ing. Using the Hartree propagators on the right-hand 
side of Eqs. ((2T|) for the self-energy we obtain a solution 
for 99 [a;+ — JI — T,{uj+)] ^ /doj = in a form 



2U{n + U\nVl + S^ )(arctan(i) -|-7Z| lna|) 
C/2 (arctan(x) + /i| lna|)2 - (tt + U\n^/l + Y 



(22) 



where we denoted x = lu /a{\ + /I^) and assumed that 
it is of order unity. In the symmetric case, Jl = 0, we 
have a trivial solution x — Q. The position of the Kondo 
resonance peak moves off the Fermi energy in the op- 
posite direction given by the effective chemical potential 
Ji = Ed + Un. For small values of Ji we explicitly obtain 



2Ufi a\ lna| 



(23) 



Maximum of the spectral function for larger effective 
chemical potentials Ji 3> 7r/|lna| eventually moves out 
from a region of order of the Kondo scale around the 
Fermi energy. 

The width of the Kondo resonance can be estimated 
from the distance between the opposite points of steepest 
descent of the spectral function around the peak near the 
Fermi energy. They can be obtained as points where the 
derivative of the real part of the one-electron propagator 
vanishes. It happens for 



{tt + U\nVTTx^ f - C/2(arctan(i) +7I| \na\ f 
2f/(arctan(i) + JI\ lna|)(7r + C/ln\/l + ) 



(24) 



For small effective chemical potentials, JI <C 7r/2|lna|, 
Eq. ([M)) has two solutions i ~ 1. The Kondo peak 
is then well formed. It broadens with increasing the 
effective chemical potential. Above the critical value 
Jl = 7r/2| lna| one of the solutions is pushed to (minus) 
infinity and the Kondo peak dissolves in the bulk of other 
states. 

The positions of the centers of the satellite Hubbard 
bands can be determined from vanishing of the imagi- 
nary part of the full one-particle propagator. We have 
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two solutions apart from the Fermi energy. The centers 
of the satellite peaks are determined from a solution pro- 
portional to the logarithm of the Kondo scale | lna|. We 
obtain explicitly for the Hartree propagators in SIAM 

_ l u\\na\ _ , ^UXee{U\Xek\~l) 

il±= H±\ — =/i±- ■. ■. (25) 

where the two-particle bubbles Xee and Xeh are deter- 
mined from Eqs. (fT51) . In the symmetric case we have 
fl± = ±{//7r. This result deviates from the atomic limit 
levels r2± — ±U/2, but we must realize that the Hub- 
bard satellite peaks in the Kondo regime are not related 
to the atomic limit, since the latter is not (cannot be) 
part of the solution with a singular two-particle vertex. 
Beware that neither the exact Bethe-ansatz solution in- 
corporates the atomic limit. Away from the half-filled 
impurity level, the distance between the satellite peaks 
increases with the square root of the effective chemical 
potential JI. Eq. ([25|l also tells us that the Hubbard satel- 
lite bands may exist only in the strong-coupling regime 
defined as J7|xe/!. | > 1- The Kondo regime sets in only 
later for f7|xe/i| ^ 1- 

VI. CONCLUSIONS 

We studied the single-impurity Anderson model in 
the strong-coupling regime with the aim to develop a 
reliable and analytically controllable approximation in 
this regime. We succeeded in simplifying the parquet 
equations to a manageable form without loosing consis- 
tency and fundamental characteristics of the Fermi-liquid 
regime for arbitrary fillings of the impurity level. We 
showed that the derived simplified parquet equations for 
two-particle irreducible vertices in the electron- hole (A'^'*) 
and electron-electron (A*^*^) singlet channels result from 
a resummation of the diagrammatic perturbation expan- 
sion in the critical region of a singularity in the two- 
particle vertex. Such a singularity emerges in SIAM due 
to multiple electron- hole scatterings. The resummation 
of the perturbation expansion proceeds so that the crit- 
ical region is reached first within RPA. Diagrammatic 
contributions beyond RPA are then summed not to be 
driven out of the critical region of the two-particle ver- 
tex. Singularity in the two-particle vertex is controlled 
by a vanishing scale. This new scale replaces the interac- 
tion strength in the perturbation expansion beyond RPA 
and is used to select relevant contributions in the critical 
regime. 

The principal assumption for rearranging the dia- 
grammatic expansion in a two-particle criticality is to 
explicitly sum and control two-particle contributions. We 
distinguish two basic classes of two-particle functions: 
singular and regular. The former diverge with vanishing 
of the critical scale and the latter remain finite. Loga- 
rithmically divergent functions (integrals over the singu- 
lar functions) are treated as (marginally) regular ones. 



The critical scale enables us to neglect finite differences 
in the regular two-particle functions and replace them 
with constants. Only the potentially divergent functions 
are kept dynamical in the two-particle criticality. Even 
more, only the singular low-frequency asymptotics mat- 
ters in the critical region. 

It is clear that contributions beyond RPA must be 
summed in a self-consistent way to make the perturba- 
tion theory of SIAM singularity-free for any finite inter- 
action strength. The parquet equations provide us with 
the necessary two-particle self-consistency. We presented 
two ways how to simplify the parquet equations in the 
critical region of the singularity in the electron-hole chan- 
nel. They differ in the way we suppress the frequency de- 
pendence of regular functions. This ambiguity, however, 
influences only non-universal critical properties. In the 
first construction we replaced the regular vertices with 
their specific averaged values. In the second one we used 
the fact that only frequencies of order of the vanishing 
critical scale remain important. In both cases we reached 
a set of manageable equations determining an effective 
interaction, replacing the bare one in RPA, and the sin- 
gular two-particle vertex. 

The solution of the simplified parquet equations led 
to the Kondo behavior in the strong-coupling regime for 
all fillings of the impurity level. We found that the 
vanishing critical scale of the two-particle vertex is the 
one determining the Kondo temperature at which the 
two-particle vertex saturates and the Kondo temperature 
overtakes the control of the low-temperature behavior of 
two-particle vertex and correlation functions. We thus 
find that the Kondo behavior is just a critical behavior 
due to a singularity in the two-particle vertex induced 
by frequency fluctuations. The critical point in SIAM is 
reached only by infinite bare interaction, but the criti- 
cal region extends deep into the physical domain of its 
finite values. Unlike the existing (approximate) solutions 
of SIAM we clearly demonstrated the two-particle origin 
of the Kondo behavior. In fact, it is not the density of 
states on the Fermi surface po that controls the Kondo 
behavior but rather the electron- hole bubble \xeh\- It is 
fortunate that for the Lorentzian DOS used here and in 
the Bethe ansatz the two numbers coincide. 

Finally we investigated one-electron functions in the 
Kondo regime. We related the two-particle Kondo scale 
with the width of the quasiparticle peak near the Fermi 
energy and found a good qualitative agreement with the 
exact result. We were able to estimate analytically the 
width of this peak and showed that it is well formed 
only close to the half-filled case. Farther away from the 
symmetric situation the Kondo resonance dissolves in the 
bulk of other states. We also showed that our parquet 
approximation not only reproduces correctly the Kondo 
resonance in the spectral function but it is also able to 
produce the satellite Hubbard bands. We found a crite- 
rion for their existence and an estimate for the positions 
of their centers. 

To summarize, we presented a simple analytic approx- 
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imation producing dynamical properties of SIAM for all 
interaction strengths and impurity fillings. It correctly 
reproduces universal features of the low-temperature 
Kondo behavior and predicts a Kondo resonance peak 
near the Fermi energy in the strong-coupling regime 
for small deviations from the electron-hole symmet- 
ric case. At intermediate couplings it is in a good 
quantitative agreement with more laborious numerical 
renormalization-group calculations. We identified the 
Kondo behavior as the critical behavior near a singular- 
ity in the electron-hole Bethe-Salpeter equation with bal- 
anced electron-hole and electron-electron multiple scat- 
terings. The electron-hole scatterings are needed to 
reach a singularity in Bethe-Salpeter equations and the 
electron-electron ones to screen the inter-particle interac- 
tion. Due to its universality, the approximation can eas- 



ily be extended to more complex multi-orbital or trans- 
lationally invariant lattice models for a reliable qualita- 
tive and quantitative investigation of a transition from 
weak to strong-coupling regimes in materials with tangi- 
ble electron correlations. 
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